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ABSTRACT 


The  method  of  difference  equation?  is  applied 
to  the  problem  of  estimating  temperature  patterns 
in  multiply-shielded  systems.  Results  are  obtained 
in  terms  of  elementary  functions  which  can  be 
readily  evaluated.  Ihe  equations  are  employed  to 
examine  the  effect  of  shield  curvature,  end 
radiation  and  gas  radiation  on  the  temperature 
recorded  by  a  shielded  thermocouple  for  a  variety 
of  operating  conditions.  Calculations  indicate 
that  whereas  shield  curvature  and  gas  radiation  may 
be  neglected  in  a  well-designed  shielded  thermo¬ 
couple,  end  radiation  can  result  in  a  substantial 
error  in  the  thermocouple- temperature  measurement. 


INTRODUCTION 


A  schematic  diagram  of  a  multiply- shielded  system  is  shown  in  Fig.  1. 

The  shield  system  is  assumed  to  be  located  in  a  flowing  gas  stream  of  temperature 
Tg  which  is,  in  turn,  contained  in  a  duct  with  a  wall  temperature  Tw. 

The  temperature  of  the  central  body  Tj_,  in  the  absence  of  any  shielding, 
may  be  computed  by  solving  a  heat-balance  equation  which  considers  radiant  heat 
exchange  with  the  duct  wall,  convection  and  radiation  to  and/or  from  the  gas  stream 
and  any  heat  generated  or  absorbed  by  the  shield.  Conduction  is  neglected  since 
this  effect  can  usually  be  made  negligible  by  suitable  insulation  of  the  supports. 

An  additional  heat  balance  equation  is  written  for  each  added  shield.  The  system  of 
equations  describing  the  shield  array  shown  in  Fig.  1  contains  Y\  unknown  temperatures 
among  the  )■}  algebraic  relations.  This  system  of  equations  may  be  solved  in  a  number 
of  ways. 

The  use  of  difference  equations  in  the  analysis  of  multiply- shielded  systems 
appears  to  have  been  neglected  in  the  engineering  literature.  This  mathematical 
procedure,  however,  yields  equations  in  terms  of  elementary  functions  which  are 
readily  evaluated.  This  caper  is,  therefore,  primarily  intended  to  present  an 
engineering  method  for  rapidly  estimating  the  temperatures  in  configurations  with 
thermal-radiation  shields. 

ASSUMPTIONS 

The  following  assumptions  are  made: 

1  The  convection  coefficient  hc  is  constant  within  the  shield  system. 

2  The  linearized  radiation  coefficient  hr  is  constant  within  the  shield 

system. 

3  The  stream  temperature  Tg  is  uniform. 

U  The  wall  temperature  Tw  of  the  duct  enclosing  the  shield  system  is 

uniform. 

5  The  temperature  of  each  shield  T^  is  uniform. 

6  Conduction  losses  are  negligible. 

The  linearized  radiation  coefficient  hr  is  generally  computed  as: 


hr  *  U  CT'  € 


where  T  is  an  "average1'  temperature  in  the  shield  system.  The  calculation  of  the 
convective  heat  transfer  coefficient  hc  and  the  "effective"  gas  stream  temperature 
Tg  has  been  amply  discussed  in  many  references  (l),l  ($)  and  is  not  included  here. 
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NOMENCLATURE 

The  following  nomenclature  is  used  in  this  paper: 

A,E,C,D  =  Constants 

E  =  end  radiation  factor 
F  »  intra-shield  radiation  factor 
G  a  gas-shield  radiation  factor 

h  *  heat  transfer  coefficient,  (Btu/hr-ft^-deg  R) 
k  *  integer  (l,  2,  ,.,n) 

Yi  a  number  of  probe  elements 
T  =  temnerature,  Deg  R. 

AT  s  Tg“T?  deg  R. 

€  -  emissivity 

CT  ~  Boltzmann  constant,  (Btu/hr-ft^-deg  R^) 

Subscripts 

c  “  convection 
g  =  gas 
i  =  innermost 
o  -  outermost 
r  a  radiation 
w  ~  wall 

1  2  o o ok 

99°  9  _  index  refers  to  shield  element 

Q  O  OU 

Other  symbols  are  defined  as  they  are  introduced  in  the  text. 
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A  SIMPLE  SHIELD  SYSTEM 


The  following  analysis  of  a  simple  shield  system  illustrates  the  general 
analytical  procedure  and,  furthermore,  leads  to  results  closely  approximating  those 
of  more  complex  systems  described  in  the  next  section- 

Consider  a  shield  system  consisting  of  a  set  of  plane  parallel  shields 
arranged  about  a  central  body  as  shown  in  Fig-  1-  Assume  that  the  external  radiation 
and  convection  heat- transfer  coefficients  are  equal  to  the  internal  radiation  and 
convection  heat* transfer  coefficients,,  respectively-  Also,  assume  that  end  radiation 
to  and  from  the  shield  assembly  is  negligible  and  there  is  no  interchange  of  radiation 
between  the  shield  system  and  the  gas  stream - 

With  reference  to  Fig-  1  the  preceding  assumptions  lead  to  the  following 
heat  balance  equation,  per  unit  shield  surface  area,  for  the  shield  element  with  index  k 


2hc(Tg  *  T}c)  +  hr(Tk-l  -  ?k)  +  M^k+l  -  Tfc)  *  0  1<  k<  n  .  .  (1) 


It  is  convenient  to  employ  the  temperature  difference  identity. 


AT  =  Tg  -  T 


in  terms  of  which  the  heat-balance  eouation  becomes 


ATi+i  -  2AATk  +  ATk-1  *0  i<k<n.  .  .  .  .(2) 

where 


Eouation  (2)  is  a  homogeneous  linear  difference  equation  of  the  second  order 
with  constant  coefficients.  The  solutions  of  difference  equations  of  this  type, 
where  A  is  a  real  constant,  are  well  known.  For  the  case  where  A>  1  the  general 
solution  is  of  the  form  (reference  2  page  2Ul<>) 


ATk  *  c1e*<k  + 


(3) 


-  h  - 


where  o C  is  defined  by  cosh  =  A,  and  Cl  and  C2  are  constants  which  are  determined  by 
the  boundary  conditions,, 

The  boundary  conditions  are  given  by  the  heat-balance  equations  written  for 
the  central  body  (k  =  l)  and  the  outermost  shield  (k  *  n)»  These  equations  are 


MTg  “  Tk)  +  MTk+l  -  Tk)  *  0  k  =  1  .  .  .  (h) 

2hc(Tg  -T^)  +  hr(Tk„i  —  Tk)  +■  hr(^k+l  **  Tk)  *  0  k  -  n  °  •  •  (5) 


In  terms  of  the  temperature  difference  A  T  the  boundary  conditions  become 

AAT^  —  A  Tj  “  0  o  •  •  o  o  (6) 

2AATn  -  A  Tn„!  »  AT,  .  .  „  .  (7) 

The  constants  C]_  and  C2  are  determined  by  substitution  of  Equation  (3)  into 
the  boundary  conditions,  Equations  (6)  and  (7),  which  yield 


G1 


At  Tw 


O 


*(8) 


Inserting  these  results  into  Equation  (3)  and  simplifying,  we  arrive  at  a 
simple  expression  for  the  temperature  of  any  element  in  the  simple  shield  system: 


ATk  „  cosh(k-l)c< 
ATW  "  cosh  n 


(10) 


The  temoerature  difference  at  the  central  body  (Tg  -  Ti)  may  be  readily 
computed  from  the  relation  obtained  by  setting  k  =  1  in  Equation  (10) 


(11) 


cosh  not 


Equation  (11)  is  convenient  for  estimating  the  number  of  radiation  shields 
necessary  to  suppress  lateral  radiant  heat  transfer  in  a  multiply-shielded  system  to 
a  desired  level  .  Eouation  (11)  is  niotted  in  Fig.  2,  where  curves  of  AT^/Z\TW  versus 
n  for  hc/hr  values  of  ?,  1  and  2  are  shown.  Mote  that  Eauation  (11)  is  valid  for  gas 
temperatures  either  greater  or  less  than  the  wall  temperature. 


THE  GENERAL  SOLUTION 

The  method  of  determining  the  temperature  pattern  in  a  shield  system  subject 
to  more  general  boundary  conditions  is  analogous  to  the  method  outlined  in  the  fore¬ 
going.  A  schematic  diagram  of  a  circular  cylindrical  multiple-shielded  system  is  shown 
in  Fig.  3.  The  heat- balance  equation  for  the  ktVl  element,  including  end  radiation 
effects,  is 


(:-'hc  *  Ohr)(Tg  =  ?k)  Fhr(Tk~l  -  Tk)  *  hr(Tk+l  -  Tk)  +  Ehr(Tw  "  Tk)  =  0 

l^fC^n  o  e  o  .  o  .  .  .(12) 


The  parameter  G  is  a  measure  of  the  amount  of  net  radiation  between  the  gas 
and  the  k^h  shield  in  relation  to  the  lateral  radiation  between  the  shields.  It 
includes  end  radiation  to  and  from  the  gas,  as  well  as  any  radiant  heat  exchange  with 
the  gas  flowing  between  the  shields. 

The  radiant- interchange  configuration  factor  F  is  the  fraction  of  radiant 
energy  leaving  the  surface  of  the  k^:  element,  which  is  intercepted  by  the  surface  of 
the  (k-l)th  element.  For  circular  cylindrical  shields 


F  c 


lk-l 


where  R  is  the  shield  radius  (h) . 

The  factor  E  is  a  measure  of  the  amount  of  end  radiation  between  the  wall 
and  the  k^h  element.  Values  for  E  may  also  be  computed  with  the  aid  of  reference  (h). 

The  parameters  E,  F,  and  G  are  assumed  constant  for  all  internal  elements 
of  the  shield  assembly.  If  they  are  considered  as  functions  of  the  index  k,  the 
coefficients  in  the  difference  eauation  are  variables,  and  no  general  solution  is 
usually  obtainable. 
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Equation  (12)  is  written  in  the  temnerature-difference  notation  a s 


2hc  +  (1  +  E  +  F  +  G)hr  m 

ATk+1 - - - r - -  ATk  *  FATk_x  =  EATW  l<k<n  .  .  .(13) 

nr 


the  are  of 


This  equation  is  reduced  to  a  "standard"  form  similar  to  Equation  (2)  by 


ATk  a  P  AUk . (1U) 


Combining  Equations  (13)  and  (lli)  we  arrive  at 


AUk+i  -  2BAUk  +  A  Vi  -  —5-  ATW  . 

„k+l 


.(15) 


where 


2hc  +  (1  +  E  +  F  +  0)hj. 

2h^ 


Eouation  (15)  is  a  nonhomogeneous  l-'near  equation  of  the  second  order  with 
constant  coefficients. 

The  general  solution  of  a  nonhomogeneous  eouation  of  the  form  of  Equation 
(15)  may  be  expressed  as  the  sum 


Auk  =  -  AtJ<p)  ....  (16) 


where  AUk  ‘  is  the  general  solution  of  the  homogeneous  eouation 


AUk+1  -  2BAUk  +  AUjf.!  =  0 


and  AkUk  is  any  particular  solution  of  Equation  (15). 
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The  general  solution  of  the  homogeneous  eruation  is 


[cf.  Eouaticn?  (2)  and 


AU^H)  =  Dxe  "k  +  h2r  ";'k  .  .  .  (17) 


where 


B  =  cosh/* 


The  narticular  solution  can  be  derived  by  the  method  of  undetermined 
coefficient?.  It  may  be  verified  that  a  narticular  solution  is 


-EAT„ 

'J - 

F2(2BF2-1-F) 


.(18) 


by  substituting  this  exoression  in  Eouation  (15). 

The  general  equation  for  A is  therefore 


AT-k  =  F 


=  ^  [Die^k 


D2  e' 


(3kJ_ 


EAT 


w 


2BF2-1-F 


1<  k  Cn 


(19) 


The  two  boundary  conditions  for  the  shield  system  shown  in  Fig.  3, 
including  end  radiation  effects,  are 


(he  +  hei  +  Gihr)(Tg  -  Tk)  +  hr(Tk+l  -  Tk)  +  Eihr(Tw  -  Tk)  =0  k  •  1  ...  (20) 


(he  *  hco  +  0ohr)(Tg  -  Tk)  +  Fhr(Tk-l  -  Tk)  +  hro(Tw  -  Tk) 

+  Eohro(Tw  -  Tk)  -  0  k  s  n . (21) 

In  the  temperature-difference  notation  the  boundary  conditions  are  written  as 

-  *T2  -  . 


(22) 


r 


JnF?ATn  -  ?ATn_i  = 


(1  +  Ec)  — -  ATW 


where 


hc  +’nci  +  (1  +  Ei  +  u})hr 


Jn  = 


hc  +  hco  +  (F  +  0o)hr  +  (1  +  E  )h 


The  constants  Dy  and  D2  in  Eouation  (19)  are  evaluated  with  the  aid  of 
the  boundary  condition  Ecuations  (22)  and  (23).  After  some  algebraic  manipulation 
the  final  eouation  for  i£ 


A  Tw 


k 


tifi  ♦  4>2'h 


lsinh(n-] 


-  2i  )sinh(n-2. 


2EF5-1-F 


=  JR  sinh(n-k)r3-  sinhCn-k-l)^ 

^2  “  J1  finh(k-l )3-  sinh(k-2)  /$ 


,  Et  ( J]_F2 -l)E 

0-,  =  —  +  - 1 - 

1  1  F  F92Err-l-F) 


<p  =  (1  +  ^)hro  t  (Jn^-  F)E 

2  F^hr  F— tj“(2EF**-1-F) 


The  refinements  of  Eouation  ( 2li )  are  inconsistent  with  the  crude  initial 
assumptions.  However,  the  eouation  may  be  used  to  determine  the  relative  imoortance 
of  several  parameters  which  have  been  assumed  to  influence  the  temperature  oattern  of 
the  shield  system. 


EFFECT  OF  ‘’FIELD  CTP,fATT.TPE  PI  TfEPrTO'TPL!:  APPLICATION'’ 

The  first  of  these  parameters  to  be  considered  is  shield  curvature.  Tie 
rati°  ^ *1  is  nlotted  as  a  function  of  the  number  of  probe  elements  in  Fig.  b  for 
several  values  of  hc/hr  with  the  configuration  factor  F  taken  as  1.0, 

0.9  and  0.8.  The  innemort  orofce  element  (k*l)  has  been  selected  for 
comparison  since  this  element  is  of  most  interest  in  shielded  thermocouple  applications 
where  curvature  effects  may  be  important. 

It  is  seen  in  Fig.  b  that  the  influence  of  sharply-curved  and/or  widely 
spaced  curved  shields  on  +he  temperature  of  the  innermost  element  of  a  shielded 
thermocouple  system  is  secondary  in  relation  to  the  influence  of  the  number  of 
shields  and  to  variations  in  convective  and  radiant  heat  transfer.  Furthermore,  the 
factor  F  is  close  to  unity  for  typical  shielded  thermocouple  geometries.  ,).  These 
assemblies  may  therefore  be  analyzed  as  parallel  shield  assembla  .-r  wit  .  on.y  higher 
order  errors  incurred  due  to  neglecting  the  effect  of  shield  c  v;  sure. 


END  RADIATION 

A  shielded  thermocouple  is  also  a  convenient  model  to  use  in  estimating 
the  contribution  of  end  radiation  to  variations  in  the  temperature  patterns  of  a 
shield  system.  The  ratio  AT^/  ATW  is  plotted  against  the  number  of  thermocouple 
elements  lq  in  Fig.  5  for  end  radiation  coefficient  (E)  values  of  0.00,  0.01  and  0.02. 
The  end  radiation  factor  E  is  of  the  order  of  0.01  to  0.02  for  typical  thermocouple 
shield  assemblies  (3). 

It  can  be  seen  in  Fig.  5  that  end  radiation  becomes  the  dominating  influence 
as  the  number  of  shields  increases.  This  is  due  to  the  fact  that  end  radiation 
depends  very  little  on  the  number  of  shields  and  becomes  completely  independent  of  the 
number  of  shields  once  lateral  radiation  is  suppressed.  For  example,  if  hc/hr  *  20 
and  E  z  0.02,  the  ratio  AT^/iiTw  is  dependent  on  the  number  of  shields,  as  follows: 


Number  of  Shields  ATl 

(n-1)  ^w 

1  .0390 

2  ,009b 

3  .00b3 

b  .003b 

5  .0033 

oo  .00327 


Fig.  5  quantitatively  illustrates  the  well-known  fact  that  end  radiation 
should  not  be  ignored  if  a  high  degree  of  accuracy  in  thermocouple  temperature  near  re- 
menta  ie  required. 
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GAf  RADIATION 

It  is  not  feasible  to  plot  the  effect  of  varying  3  for  a  thermocouple 
since  the  curves  for  various  values  of  G  tend  to  overlap.  this  effect  i:  illustrated 
ty  the  following  tabulated  results: 


&  T]_ 


Number  of  fhields 

Value 

of  G 

(n-1) 

0.00 

0.01 

1 

.091:8 

.0911 

2 

.0271 

.0269 

3 

.0089 

.0089 

h 

.  00li06 

.001:06 

Conditions  -  hc  = 

hpi  u  hco  8 

hr  =  hro 

E  —  Ej[_  *  E0  -  0.01 
F  =  1.0 


As  the  number  of  shields  increases  the  temperature  of  the  probe  approaches 
the  gas  temoerature  and  no  appreciable  radiant  energy  is  exchanged  between  the 
thermocouple  and  the  gas  stream. 


CONCLUSIONS 

The  preceding  analysis  presents  an  engineering  method  for  rapidly  estimating 
the  temperatures  in  a  system  which  contains  an  array  of  thermal  radiation  shields. 

The  application  of  the  calculus  of  finite  differences  leads  to  an  equation  which  may 
be  used  to  compute  the  temperature  of  any  shield  once  the  system  parameters  have  been 
evaluated.  The  method  applied  to  the  problem  of  tempera' ure  measurement  by  shielded 
thermocuples  produces  results  which  are  Qualitatively  well  known,  e.g.  shield  curvature 
and  gas  radiation  effects  are  secondary  quantities  in  a  well- shielded  thermocouple; 
end  radiation  ma-  be  important.  The  Quantitative  influence  of  convection,  radiation, 
the  number  of  shields,  etc.  on  the  thermocouple  reading  mav  be  seen  readily  from  the 
graphs  which  are  included. 
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